Energy  Conversion  and  Management  87  (2014)  495-503 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Energy  Conversion  and  Management 

journal  homepage:  www.elsevier.com/locate/enconman 


Optimization  of  a  thermoelectric  system  for  power  generation  realized 
by  “hidden”  components 


CrossMark 


Maria  La  Gennusa,  Massimo  Morale,  Gianluca  Scaccianoce  * 

Dipartimento  di  Energia,  Ingegneria  dell’Informazione  e  Modelli  Matematici,  Universita  degli  Studi  di  Palermo,  Viale  delle  Scienze  9,  Palermo,  Italy 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  21  November  2013 
Accepted  10  July  2014 
Available  online  3  August  2014 


Keywords: 

Thermoelectric  effect 
Solar  thermal  collector 
Peltier  cell 


Today  the  utilization  of  renewable  energy  sources  is  very  important  to  support  the  traditional  ones,  and 
also  in  order  to  reach  the  objectives  fixed  by  the  Kyoto  Protocol. 

In  this  context,  also  well-known  technologies  related  to  the  electric  energy  production,  formerly  not 
suitable  -  due  to  the  poor  technological  level  of  the  components  or  due  to  the  high  costs  linked  to  the 
system  efficiency  -  may  now  find  a  suitable  collocation. 

Among  these  technologies,  the  thermoelectric  effect  may  be  very  advantageous  in  some  particular 
applications  of  small  power.  A  power  generation  system  utilizing  this  technology  has  been  studied  by 
the  authors.  This  system  utilizes  the  thermal  difference  between  a  lower  temperature  sink,  e.g.  a  water 
reservoir,  and  a  higher  temperature  source,  such  as  exhaust  heat  (e.g.  from  industrial  processes)  or  as 
produced  by  solar  energy  captured  by  an  extended  surface.  The  thermoelectric  generator  is  made  up 
of  thermocouples  placed  on  the  internal  surface  of  a  tubular  heat  exchanger,  located  inside  the  water 
body  and  filled  up  with  hot  water  coming  from  the  energy  source. 

In  this  paper,  the  physical  mechanisms  and  the  energy  balances  governing  the  system  are  described  in 
order  to  evaluate  its  energetic  feasibility.  Once  fixed,  the  electrical  power  required  by  the  end  user,  the 
system  parameters  are  optimized  with  the  aim  to  have  the  most  compact  equipment  and  components. 
A  simulation  of  the  proposed  system  was  performed  referring  to  the  city  of  Palermo,  Italy.  A  very  poor 
economic  convenience  of  the  system  has  been  found  in  this  applications. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

A  greater  sensitiveness  in  public  opinion  towards  environmen¬ 
tal  problems,  highlighted  by  the  media,  especially  in  the  more 
economically  developed  countries,  joined  with  an  increasing  mean 
cost  of  fossil  fuels,  results  in  a  more  intensive  utilization  of  renew¬ 
able  energy.  This  is  also  to  reach  the  objectives  signed  by  countries 
in  the  Kyoto  protocol  [28]. 

In  the  wide  range  of  renewable  energies  in  order  to  supply 
electrical  energy,  Research  and  Development  look,  with  a  new 
interest  at  techniques  previously  well-known  but  practically  not 
implemented  due  to  their  low  efficiencies,  related  either  to  the 
technological  level  of  components,  or  to  the  unreasonable  costs. 
Among  these  techniques  there  is  the  thermoelectric  effect,  charac¬ 
terized  by  low  efficiencies  but,  in  some  particular  applications, 
suitable  to  give  good  advantages:  e.g.  in  small  electric  power  appli¬ 
cations  where  a  regular  electric  supply  is  not  suitable. 
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The  thermoelectric  elements,  or  thermocouples,  exploit  the 
principle  observed  by  the  German  scientist  Thomas  Johann 
Seebeck  in  1821  and,  separately,  observed  by  the  French  scientist 
Jean  Charles  Athanase  Peltier  in  1834.  In  a  circuit  made  of  two 
different  metals  by  placing  the  junctions  at  two  different  temper¬ 
atures,  an  electromotive  force  is  generated.  Vice  versa  by  applying 
an  electromotive  force,  it  is  possible  to  have  different  temperatures 
at  the  junctions  and  so  a  heat  exchange.  This  effect  is  widely  used 
in  many  engineering  applications;  it  is  the  basic  principle  for  a  very 
easy  to  use,  reliable  and  versatile  temperature  measuring  system 
[23]. 

The  system  analyzed  in  this  study  [16]  uses  the  thermal  differ¬ 
ence  between  a  reservoir  at  a  low  temperature  (e.g.  groundwater, 
river),  and  a  reservoir  at  a  high  temperature  (e.g.  waste  heat  in 
an  industrial  process,  heat  captured  by  a  solar  collector). 

The  system  was  previously  designed  to  be  used  as  an  aid  for  the 
development  of  African  villages,  where  an  electrical  supply  was  not 
available.  The  system  components  may  be  considered  “hidden” 
because  the  heat  exchanger  at  a  high  temperature  may  be  placed 
on  a  roof  (solar  collector)  or  under  the  ground.  The  other  heat 
exchanger  at  a  low  temperature  must  be  dipped  in  a  water 
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Nomenclature 

Latin  symbols 

P 

volumetric  thermal  expansion  coefficient,  K-1 

Ac 

overall  surface  area  of  solar  collectors,  m2 

S 

density,  kg/m3 

Chw 

Hazen-Williams  coefficient 

S * 

solar  declination,  rad 

cp 

specific  heat  at  constant  pressure,  J  kg-1 1<-1 

&pn 

difference  between  the  Seebeck  coefficients  of  p-  and 

D 

diameter,  m 

n-type  elements,  V/K 

d 

thickness,  m 

4> 

latitude,  rad 

dglass 

thickness  of  single  cover 

7 

tilt  angle,  rad 

E 

ends  potential  difference,  V 

de 

pump  electrical  efficiency 

f,fr 

solar  collector  factors 

Pfin 

fin  efficiency 

h 

convective  heat  transfer  coefficient,  W  m-2  K-1 

Pm 

pump  mechanical  efficiency 

I 

irradiance,  W  m-2 

cp ,  6 

angle,  rad 

Kglass 

extinction  coefficient  of  glass,  m_1 

2 

thermal  conductivity,  W/(m  K) 

/,  L 

length,  m 

li 

viscosity,  Pa  s 

M 

mass  flow  rate,  kg/s 

V 

kinematic  viscosity,  m2/s 

m 

solar  collector  useful  parameter 

P 

resistivity,  Q  m 

N 

number  of  cover  of  solar  collector 

Pg 

albedo 

n 

thermocouple  density  or  cover  refractive  index 

T 

transmittance  coefficient  of  glass 

P 

power,  W 

CD 

hour  angle,  rad 

Pw 

pump  power,  W 

c 

minor  pressure  loss  coefficient 

p 

pressure,  Pa 

(t«) 

transmittance-absorptance  product 

R 

electrical  resistance,  Q 

Rb 

ratio  of  beam  irradiance  on  the  tilted  surface  to  that  on 

Sub  and  superscripts 

the  horizontal  surface 

cp 

circumferential 

r 

unpolarised  radiation 

b 

beam 

S 

solar  radiation  absorbed  by  the  solar  collector,  W  m-2 

c 

cold  fluid 

T 

temperature,  °C 

d 

diffuse 

Qu 

heat  flux,  W 

e 

external 

u 

velocity,  m/s 

h 

hot  fluid 

u,uL 

thermal  transmittance,  W  m-2  K_1 

i 

internal 

w 

step  between  tubes  in  the  solar  collector,  m 

n 

n-  type  thermoelement  or  normal  direction 

X 

linear  abscissa,  m 

out 

outdoor 

z 

figure  of  merit,  K_1 

P 

p- type  thermoelement 

ZT 

dimensionless  figure  of  merit 

s 

surface 

sc 

solar  collector 

Greek  symbols 

X 

axis  direction 

ot 

absorbance  coefficient  of  slab  of  solar  collector 

reservoir.  Therefore  they  are  difficult  to  reach,  so  obtaining  the  goal 
to  protect  them  from  vandalism  and  theft. 

In  this  work,  starting  from  this  design  hypothesis,  the  system  is 
analyzed  with  the  aim  to  highlight  the  most  relevant  design 
parameters. 

Riffat  and  Xiaoli  Ma  [21]  showed  the  convenience  of  thermo¬ 
electric  effect  in  order  to  produce  electric  power  and  note  that  a 
thermoelectric  system  is  adequate  in  cases  where  thermal  energy 
is  available  at  a  low  cost  or  free.  On  the  other  hand,  they  also  note 
that  these  systems  have  a  relatively  low  efficiency,  typically 
around  5%.  Currently  many  researchers  are  working  on  thermo¬ 
electric  generators  mainly  referring  to  potential  applications 
[18,20,30,11]  and  mathematical  and  physical  characterization 
[15,10,22].  Furthermore,  several  research  works  have  been  focused 
on  using  solar-thermoelectric  systems,  also  combined  with  photo¬ 
voltaic  or  another  technologies  [17,4,31,13,29,5]. 

Many  works  in  literature  show  how  R&D  focus  on  new  materi¬ 
als  with  a  higher  efficiency  and  with  a  larger  figure  of  merit 
[19,26,7,8,24,12]. 

2.  The  proposed  model  system 

The  proposed  system  is  basically  constituted  of  four  sections: 
(1)  the  first  section  is  composed  of  the  electric  user  system; 


(2)  the  second  section  is  composed  of  the  hydraulic  circuit;  (3) 
the  third  section  is  composed  of  the  heat  exchanger  of  the  power 
source  (for  example,  solar  collectors);  (4)  the  fourth  and  principal 
section  referring  to  the  heat  exchanger  dipped  in  a  cold  reservoir 
(the  heat  sink),  that  is  also  the  thermoelectric  generator.  Fig.  1 
shows  a  simplified  sketch.  The  electrical  power  produced  by  ther¬ 
moelectric  generator  has  to  meet  the  requirements  of  the  electrical 
load  of  the  considered  site. 

For  convenience,  we  directly  report  the  main  general  relation¬ 
ships  used  in  our  simulation  model. 

The  fluid  speed  in  the  circuit  can  be  calculated  by  means  of  the 
mass  flow  rate: 

U  =  4M/(d%D])  (1) 

where  M  is  the  flow  rate,  S  is  the  fluid  density  and  Dz  is  the  inner 
diameter. 

The  major  losses  are  computed  with  the  Hazen-Williams  equa¬ 
tion  [2]: 

/  u  Y 1852  /1Y 1167 

APm,or  =  6.8l9L(— )  (-)  (*■*)  (2) 

where  CHW  is  the  coefficient  of  roughness,  L  is  the  length  of  the  pipe, 
and  g  is  the  earth’s  standard  acceleration  due  to  gravity. 
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Fig.  1.  A  sketch  of  the  proposed  system. 


The  minor  losses  can  be  evaluated  with  reference  to  the  partic¬ 
ular  coefficient  C  of  the  components  that  causes  the  pressure  loss, 
i.e., 


The  thermal  diffusivity  and  the  kinematic  viscosity  were 
calculated  utilizing  their  definition  as  the  function  of  the  other 
computed  variables. 


A V 


minor 


(3) 


For  the  evaluation  of  the  internal  heat  transfer  coefficient  of  the 
tube,  hi,  it  is  possible  to  use  the  well-known  Dittus-Boelter  equa¬ 
tion  14],  with  the  limits:  0.7  <  Pr  ^  160;  ReD  >  10,000;  D/L  >10: 

Nud  =  0.023Re;f  Pr2  (4) 


where  z  =  0.4  for  heating  (Ts  >  Tfiuid )  and  z  =  0.3  for  cooling 
(Ts  <  Tfiuid),  being  Ts  and  Tfiuid  temperatures  of  the  surface  of  the  tube 
and  of  the  fluid  that  flow  through,  respectively. 

For  the  evaluation  of  the  external  heat  transfer  coefficient  of  the 
tube,  he,  it  is  possible  to  use  the  Churchill-Chu  correlation  [9],  with 
the  limit  10-5  <  RaD  ^  1012: 


NuD  =  0.60  + 


_ r 

[1  +(0.559/Pr)9/16]8/27J 


0.387Ray6 


(5) 


In  order  to  evaluate  the  dimensionless  numbers,  we  have  to 
compute  the  thermo-physical  characteristics  of  the  fluids  (water) 
for  the  given  value  of  the  temperature. 

The  thermal  conductivity  2  and  the  density  3  of  the  water  may 
be  evaluated  utilizing  the  equations  derived  by  Choi  and  Okos  [1], 
under  the  limit  -40  °C  ^  T^  150  °C: 


2  A.  The  electrical  load 

The  electrical  load  can  be  a  generic  system,  for  example  a  water 
pump  for  a  well.  Obviously  the  design  of  the  proposed  system 
depends  on  the  electrical  power  needed  by  the  user  and  by  the 
available  heat  source.  It  is  worth  noting  that,  apart  from  the 
start-up  of  the  system,  it  is  a  stand-alone  one,  it  is  possible  to  com¬ 
plete  the  plant  with  an  array  of  batteries  in  order  to  compensate 
the  lack  of  a  power  source. 

2.2.  The  hydraulic  circuit 

A  length  of  well  insulated,  diameter  fixed  copper  pipe  and  a  cir¬ 
culation  pump  constitute  the  hydraulic  circuit.  The  pump  sends  the 
hot  fluid  through  the  pipe  and  between  the  two  heat  exchangers. 
As  parameters,  we  also  have  the  electric  power  of  the  pump  and 
the  pressure  loss  along  the  hydraulic  circuit. 

After  setting  up  the  geometry  of  the  plant,  it  is  possible  to  eval¬ 
uate  the  mass  flow  rate  of  the  hot  fluid  and  then,  by  means  of  the 
Bernoulli  equation,  to  compute  the  pressure  losses. 

The  pump  power  Pw  is  related  to  the  mass  flow  rate  of  the  hot 
fluid,  pressure  losses  and  moreover  with  mechanical  and  electrical 
efficiencies  of  the  pump.  So  we  have: 


X  =  5.7109  ■  KT1  +  1.7625  •  10“3I  -  6.7036  ■  10“6T2  (6) 

§  =  9.9718  •  102  +  3.1439  •  10“3I  -  3.7574  •  10"3r2  (7) 

The  specific  heat  cp,  the  viscosity  p  and  the  expansion  coeffi¬ 
cient  p  of  the  water  are  evaluated  by  means  of  the  following 
equations,  which  we  obtained  by  means  of  an  interpolation  of 
the  related  data  tables  9],  under  limit  0  °C  ^  T  ^  120  °C: 

cp  =4207.7-  1.06181  + 0.01 14T2  (8) 

p  =  1.738  •  1CT3  -  5.020  ■  10"5r  +  7.961  ■  KT7!2  +  6.399 

•KT9T3  + 1.992  KT11!4  (9) 

P  =  -6.129  •  1CT5  +  1.613  ■  KT5T  -  1.594  ■  W7T2  +  1.042 

.  10"9r3- 2.501  10  12T4  (10) 


Pw  =  MAp/(St]mt]e)  (11) 

where  A p  indicates  the  pressure  losses  in  the  whole  system,  the 
sum  of  major  and  minor  losses:  Apmajor  +  Apminor,  as  stated  above. 

2.3.  The  heat  exchanger  of  the  heat  source 

In  this  work  we  assume  that  the  heat  exchanger  of  the  heat 
source  is  a  solar  collector  and  the  main  output  parameters  are 
the  pressure  loss  of  the  solar  collector  and  outlet  temperature  of 
the  fluid  for  fixed  flow  rate  and  climatic  parameters. 

2.3 A.  The  mathematical  model  of  the  solar  collector 

Referring  to  the  half  width  of  the  sheet  (the  length  of  fin)  for 
each  tube  is  calculated  by  means  of  the  following  equation: 


2 


(12) 
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where  W  is  the  width  of  the  step  between  tubes  and  Dsc>e  is  the 
outer  diameter  of  the  tube  in  the  collector. 

It  is  now  useful  to  introduce  a  parameter  m  that  can  be  utilized 
in  the  following  relationship: 


The  pertinent  incident  angle,  9 ,  is  different  in  the  case  of  beam, 
diffuse  or  reflected  radiation.  Duffie  and  Beckman  [6]  provide  the 
following  equation  to  calculate  these  incident  angles: 

9b  =  arccos[cos (</>  -  y)  cos  <5*  cos  co  +  sin(</>  -  y)  sin  3*]  (23) 

9d  =  59.7  -  0.1388y  +  0.001497y2  (24) 


where  UL  is  the  thermal  transmittance  of  the  solar  collector,  A  is  the 
thermal  conductivity  of  the  slab  and  d  is  the  thickness  of  the  slab. 

The  function  of  the  standard  fin  efficiency  for  straight  fins  with 
a  rectangular  profile  is: 

tanh (mj) 

m.|  (14) 


The  collector  efficiency  factor  is  provided  by  the  following 
equation: 


F'  = 


W- 


Dsc,e+2lrjjjn  7lDf 


:DSc,hj  J 


(15) 


Aground  =  90.0  -  0.5788y  +  0.002693r2  (25) 

Instead,  the  Rb  coefficient  is  provided  by  the  following  equation: 

_  cos(</>  —  y)  cos  3*  cos  oo  +  sin(</>  —  y)  sin  3* 
b~  cos(</>)cos<f  cosco  +  sin(</>)sin<f 

where  4>  is  the  latitude  of  site,  3 *  is  the  solar  declination,  co  is  the 
hour  angle. 

In  conclusion,  the  outlet  temperature  is  provided  by: 


T  outlet  —  T  inlet  + 


Qu 

M-Cn 


(27) 


where  DSCii  is  the  inner  diameter  of  the  tube  of  the  solar  collector 
and  hi  is  the  internal  convective  heat  transfer  coefficient  provided 
by  Eq.  (4),  where  n  has  been  fixed  at  the  value  of  0.4. 

The  removal  factor  of  the  solar  collector  is  given  by: 


while  the  pressure  drop  caused  by  fluid  friction  is  provided  by 
Eqs.  (2)  and  (3). 

2.4.  The  thermoelectric  generator 


Fr  = 


M  •  Cp 

AUI 


(16) 


where  M  is  the  mass  flow  rate,  cp  is  the  specific  heat  at  a  constant 
pressure  and  Ac  is  the  overall  surface  area  of  solar  collectors. 

The  following  equation  provides  the  amount  of  the  actual  useful 
energy  gain  [6] 


Qu  =AcFR[S(m)n  -  UL(Ji  -  Tout)]-  (17) 

where  S  is  the  solar  radiation  absorbed  by  the  solar  collector,  (toc)„ 
is  the  transmittance-absorptance  product  for  normal  direction,  Tt  is 
the  inlet  temperature  of  the  hot  fluid  and  Tout  is  the  outside  air  tem¬ 
perature.  The  parameter  5  is  obtained  by  the  following  equation: 

S  =  IbRb( m)b  +  Id( m)A  +  Pg(/b  +  Id)  {'£&) ground  ~ ^ 

(18) 


where  Ib  and  Id  are  beam  and  diffuse  components  of  irradiance,  Rb  is 
the  ratio  of  beam  on  the  tilted  surface  to  that  on  the  horizontal  sur¬ 
face,  y  is  the  tilt  angle,  pg  is  the  diffuse  reflectance  of  the  surround¬ 
ings  (albedo)  and  T0tb ,  T0td  and  Tocground  are  the  transmittance- 
absorptance  products  for  different  incident  angles  that  are  provided 
by  Duffie  and  Beckman  [6].  They  assume  the  generic  form: 

(TOC)  =  1.01  T  OC  (19) 

where  t  is  calculated  by  the  following  equations: 


t  =  exp 


N  •  K glass  '  dgiass\  1/  1  ^  1  ?  \  \ 

cos^  J  '2  Vl  +  (2 N  -  l)r,|  1  +  (2 N  -  1  )r  J 

(20) 


where  02  is  the  angle  of  refraction  (arcsin()  of  ratio  of  sin()  of 
pertinent  incidence  angle,  6 ,  to  refractive  index,  n),  N  is  the  number 
of  glass  covers,  Kgiass  is  the  extinction  coefficient,  dgiass  is  the  thick¬ 
ness  of  the  single  cover,  r\\  and  r±  are  parallel  and  perpendicular 
components  of  the  unpolarised  radiation  respectively: 


r±  = 


r  M  = 


sin2  (02  -  9) 
sin2  (02  +  9) 

(21) 

tan2  (02  -6) 
tan2  (02  +  9) 

(22) 

The  model  of  the  thermoelectric  generator  was  fully  developed 
by  Suzuki  and  Tanaka  [27],  referring  to  the  case  of  a  single  tube  and 
isothermal  conditions  of  the  cold  fluid. 

The  thermoelectric  generator  is  the  main  section  of  the  pro¬ 
posed  system;  it  is  made  up  of  a  heat  exchanger  with  a  set  of  ther¬ 
mocouples,  dipped  in  the  cold  reservoir  (the  heat  sink).  Basically,  it 
is  represented  by  a  cylindrical  tube  coated  with  a  single  layer  of 
thermocouples  (Fig.  2).  It  has  also  been  assumed  that  the  p-  and 
n-type  semiconductors  are  homogeneously  aligned  perpendicular 
to  the  heat  flow,  without  an  open  space,  and  connected  electrically 
in  series.  This  is  an  ideal  configuration  of  the  thermoelectric 
generator. 

The  overall  heat  transfer  coefficient  U  of  the  heat  exchanger 
dipping  in  the  reservoir,  considering  one  unit  of  tube  length,  is 
given  by: 


U  = 


2n 


1  ,  ln(De/Dj)  ,  1 

hiPi/ 2  _r  2  _r  heDe/2 


(28) 


where  A  is  the  thermal  conductivity  of  thermocouples  defined  by: 


Ap(pp  +  Ancpn 
(pp  +  Vn 


(29) 


Fig.  2.  A  rough  sketch  of  a  section  of  the  thermoelectric  generator. 
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In  the  previous  equation  Xp  and  2n  are  the  thermal  conductivity 
of  each  type  of  junction,  and  (pp,  cpn  are  the  angles  respectively  of 
p-  and  n-type  semiconductors  of  thermocouples  (Fig.  2). 

Fixing  a  reference  x  on  the  axis  of  the  tube  and  putting  the  0 
position  at  the  beginning  of  the  tube,  the  temperature  trend  of 
the  hot  fluid  along  the  tube  is  given  by: 

T(x)  =  (Th  -  Tc)e('^)  +  Tc  (30) 

where  Th  is  the  temperature  of  the  hot  fluid  at  the  inlet  section  and 
Tc  is  the  temperature  of  the  water  in  the  reservoir. 

In  order  to  evaluate  the  electrical  potential  difference,  it  is  nec¬ 
essary  to  know  the  surface  temperature  of  the  hot  and  cold  sides  of 
the  thermocouples,  given  by: 

W  -  ™  -  5®  <T«  -  ™  T«<*>  -  ■ r‘  ■ +  as 1 nx) - ' r‘> 

(31) 


Therefore  it  is  possible  to  calculate  the  electrical  potential 
difference: 


E  = 


tl(ptlx£pn 


rs,c(x)]dx 


(32) 


where  and  nx  are  respectively  the  number  of  thermoelectric 
pairs  along  the  circumference  and  the  number  of  thermoelectric 
pairs  along  the  axis  direction,  L  is  the  length  of  the  thermoelectric 
generator. 

It  is  necessary  to  consider  the  internal  electrical  resistance 
given  by: 


R,  =  n^L 


(33) 


The  output  electrical  power,  optimized  by  balancing  the  inter¬ 
nal  and  external  resistances,  is  given  by: 


p  _  _  n (p nx [Mcp SpnjTh  Tc)]  (In (De/Dj)  r  _ 

~4Ri~  16(RP  +  Rn)L  [  nA  l 


where  Rp  and  Rn  are  respectively  the  resistance  of  p-  and  n-type 
elements. 


(35) 


The  end  of  maximizing  the  output  power  has  been  achieved  by 
adopting  the  following  angles  of  p-  and  n-type  elements: 


(36) 


3.  Simulation 

With  the  aim  to  evaluate  the  performances  of  the  proposed  sys¬ 
tem,  the  authors  have  built  a  simulating  software,  in  C-language, 
that  utilizes  the  above  mathematical  model.  Some  input  values 
regarding  the  thermoelectric  generator  were  taken  from  Suzuki 
and  Tanaka  [27].  In  order  to  assess  the  electrical  yearly  energy  pro¬ 
duction,  meteorological  data  of  a  whole  year  were  utilized  as  input 
data  in  order  to  evaluate  the  thermal  production  of  the  solar 
collectors. 

The  simulation  of  the  performances  of  the  system  has  been  con¬ 
ducted  in  a  quasi-steady  state  regime.  This  choice  directly  depends 
on  the  structure  of  the  proposed  model  that  describes  the  thermo¬ 
physical  behaviour  of  the  system  in  a  time-invariant  way.  This 
allows  to  lead  to  a  numeric  step  by  step  simulation  where,  in  each 
analyzed  step,  the  system  should  reach  its  regime  conditions.  Due 


to  the  structure  of  the  working  cycle  and  depending  on  the  fluid 
circulating  through  it,  the  time  step  selected  for  the  simulation 
was  an  hour.  The  kind  of  climatic  data  set  utilized  as  an  input  for 
the  simulation  is  the  Test  Referee  Years  (TRY),  which  provides 
hourly  values  of  climatic  parameters  of  the  selected  sites.  If  TRY 
is  not  available,  it  is  possible  to  use  a  reduced  weather  dataset  [3]. 

The  utilized  climatic  database  is  the  Test  Reference  Year  for  the 
city  of  Palermo,  Italy  [25];  Figs.  3  and  4  show  daily  trends  of  air 
temperature  and  solar  radiation  respectively.  Other  input  values 
were  taken  from  technical  and  commercial  literatures. 

3.1.  Values  of  parameters  used  in  the  simulation 

It  was  assumed  that  the  collectors  were  installed  in  the  city  of 
Palermo  (38.12N,  13.34E  -  South  of  Italy),  oriented  South,  with  a 
tilt  angle  equal  to  30°,  and  the  albedo  of  the  surrounding  area 
was  equal  to  0.2. 

The  main  characteristics  of  the  solar  collector  are  reported  in 
Table  1,  while  Tables  2  and  3  report  the  main  parameters  of 
thermoelectric  elements  and  thermoelectric  generator  respec¬ 
tively.  The  figure  of  merit  of  the  thermoelectric  module,  Z,  is 
2.59  10-3  K_1. 

The  number  of  thermocouples,  n^,  in  a  circular  section  is  calcu¬ 
lated  by  considering  a  constraining  of  the  minimum  width  of  6  mm 
for  each  thermocouple. 

The  main  parameters  of  the  thermo-hydraulic  circuit  used  are 
reported  in  Table  4. 

4.  Results  and  discussion 

After  setting  the  above  system  parameters,  running  the  simulat¬ 
ing  software  provides  the  yearly  data  output.  The  software  gives  as 
output  the  net  electrical  energy  production,  the  efficiency  of  the 
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Fig.  3.  The  trend  of  max  and  min  daily  air  temperature  of  Palermo’s  TRY. 
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Table  1 

Main  characterizing  values  of  the  solar  collector. 


Parameter  Value 


Surface,  Ac  ( variable  parameter ) 

2-20  m2 

Overall  loss  coefficient,  UL 

8.0  W  m-2  K-1 

Number  of  glass  covers,  N 

1 

Thickness  of  glass,  dg/ass 

2.3  mm 

Extinction  coefficient,  Kgiass 

16.1 

Cover  refractive  index,  n 

1.526 

Transmittance-absorptance  product,  (toc) 

0.826  584 

Riser  tube  spacing,  W 

0.15  m 

Riser  tube  number 

5 

Length  of  each  tube 

2  m 

Outer  diameter  of  the  riser  tube,  Dsc  e 

10  mm 

Inner  diameter  of  the  riser  tube,  Dsc>i 

9.6  mm 

Absorber  sheet  thickness,  d 

0.5  mm 

Thermal  conductivity,  A 

385  W  m-1  K1 

Table  2 

Characterizing  values  of  the  thermoelectric  elements. 

Thermoelectric  element  Absolute  Seebeck 

Electric 

Thermal 

coefficient  s 

resistivity  p 

conductivity  A 

(HVK-1) 

(pQ  m) 

(W  ITT1  K”1) 

p-Type  element  162 

Bi-54.3  at.%  Te 

5.55 

2.06 

n-Type  element  -240 

Bi-64.5  at.%  Te 

10.1 

2.02 

Table  3 

Characterizing  values  of  the  thermoelectric  generator. 

Section  Parameter 

Value 

Thermoelectric 

Length,  L 

10m 

device 

Thickness  of 

0.02  m 

thermocouples 

Number  density  of 

50  ITT1 

thermocouples,  nx 

Inner  diameter,  D, 

9.6,  14.0,  18.0,  23.2,  29.0, 

( variable  parameter ) 

36.20,  45.80  mm 

Cold  fluid 

Temperature,  Tc 

15  °C 

(Water) 

Table  4 

Characterizing  values  of  the  thermo-hydraulic  circuit. 


Section 

Parameter 

Value 

Pipe 

Type  of  tube 

DN90,  PE80-PFA8 

Length  of  pipe 

40  m 

Coefficient  of  roughness,  CHW 

150 

Inner  diameter,  Di  ( variable 

9.6,  14.0,  18.0,  23.2,  29.0, 

parameter) 

36.20,  45.80  mm 

Working 

fluid 

Speed,  u  ( variable  parameter) 

0.3-2.3  m/s 

Pump 

Electrical  efficiency,  rje 

95% 

Mechanical  efficiency,  rjm 

80% 

whole  system,  the  heat  exchanged  by  the  thermoelectric  generator 
and  the  sink,  and  the  electrical  power  produced.  This  provides  the 
average,  minimum  and  maximum  values  of  the  above  parameters. 
The  net  electrical  power  is  defined  as  the  electrical  power  produc¬ 
tion  (see  Eq.  (34))  minus  the  pumping  power  consumption  (see  Eq. 

(ID). 

Each  of  these  values  was  computed  varying  the  working  fluid 
velocity,  the  pipe  diameter  (both  thermoelectric  generator  and 
hydraulic  circuit)  and  the  surface  of  the  solar  collector  system. 

Fig.  5  shows  that,  in  the  range  of  the  speed  considered  in  the 
simulations  (0.3-2.3  m/s),  the  trend  of  the  net  electrical  production 
has  a  maximum  for  each  diameter.  Moreover,  the  number  of  curves 


that  show  a  maximum  value  seems  to  increase  with  the  increasing 
of  the  collector  surface,  and  the  maximum  occurs  for  the  lower 
speed  of  the  working  fluid  when  the  diameter  increases.  Further¬ 
more,  the  optimal  value  of  the  working  fluid  speed  is  around 
0.8  m/s.  In  the  case  of  high  diameters,  the  system  is  oversized 
and  requires  a  larger  collector  surface,  so  hereinafter  we  refer  only 
to  the  9.6  mm  diameter,  because  it  appears  to  be  in  compliance 
with  the  aim  of  the  application  of  the  case  study. 

The  plots  of  the  Net  Energy  Output  varying  the  speed  of  work¬ 
ing  fluid  are  reported  in  Fig.  6,  for  selected  solar  collector  surfaces. 
As  expected,  the  diagram  shows  that,  increasing  the  solar  collector 
surface,  Net  Output  Energy  grows.  It  also  highlights  that  there  is  a 
maximum  for  each  curve  and  that,  increasing  the  solar  collector 
surface,  the  maximum  value  occurs  for  a  slightly  higher  speed. 

Similarly  in  Fig.  7  the  maximum  efficiency  of  the  whole  system 
has  the  same  behaviour.  The  more  the  solar  collector  surface 
increases,  the  more  the  maximum  efficiency  of  the  whole  system 
grows.  The  increment  of  the  system’s  maximum  efficiency 
decreases  with  the  growing  of  the  solar  collector  surface  (lines 
are  closer  increasing  collector  surface).  For  a  lower  value  of  the 
working  fluid  speed,  Fig.  7  also  shows  that  it  is  possible  to  achieve 
a  better  efficiency  with  a  smaller  solar  collector  surface. 

Fig.  8  shows  that  the  maximum  electrical  power  of  the  system 
has  similar  behaviour  to  that  of  the  maximum  efficiency  with  the 
difference  that,  increasing  the  speed  values,  the  slope  of  the  curves 
decreases. 

The  maximum  temperature  drop  of  the  working  fluid  decreases 
exponentially  with  the  increasing  of  the  working  fluid  speed  and, 
considering  the  corresponding  increasing  of  the  solar  collector 
surface,  the  starting  slope  of  curve  is  more  accentuated,  as  shown 

in  Fig.  9. 

In  Fig.  10  the  trend  of  the  daily  mean  temperature  difference 
between  the  two  sides  of  the  TE  generator  are  reported  in  the  year 
for  three  different  speeds,  for  D*  =  9.6  mm  and  Ac  =  8  m2.  These 
plots,  especially  for  the  lower  speeds,  reflects  the  trend  of  daily  glo¬ 
bal  solar  radiation:  values  oscillate  with  an  amplitude  that 
amounts  to  15  °C  for  the  lowest  speed  of  0.3  m/s,  to  a  value  of 
5  °C  for  the  highest  speed  of  2.3. 

It  is  also  possible  to  see  that,  increasing  the  speed  value  of  the 
circulating  fluid,  curves  run  with  a  higher  value  of  temperature. 
From  a  mean  AT  of  around  12.5-20  °C  for  u  =  0.3,  in  the  case  of 
u  =  2.3  the  mean  value  is  around  40  °C.  Moreover  the  dimension¬ 
less  figure  of  merit,  ZT,  varies  from  0.750  to  0.826. 

In  Fig.  1 1  the  temperature  difference  between  the  two  sides  of 
the  TE  module  versus  hourly  outside  air  temperature  are  reported, 
for  Di  =  9.6  mm  and  Ac  =  8  m2.  Basically  it  is  possible  to  see  that  the 
values  of  outside  air  temperature  are  concentrated  in  the  range  of 
10-30  °C  and  that  there  is  no  particular  dependence  between  the 
two  parameters.  It  is  noteworthy  that  increasing  speed,  as  stated 
in  Fig.  1  la-c,  the  minimum  values  of  ATcoid-hot  occur  for  the  higher 
temperature.  This  is  due  to  higher  losses  that  do  not  permit  the 
imposed  self-reliance  of  the  system.  Furthermore,  the  higher 
values  of  ATcoid-hot  occur  for  the  higher  speed  of  working  fluid. 

Fig.  12  shows  that  the  temperature  difference  between  the  two 
sides  of  the  TE  modules  depend  on  global  solar  radiation,  for 
Di  =  9.6  mm  and  Ac  =  8  m2.  The  more  global  solar  radiation 
increases,  the  more  the  temperature  difference  between  the  two 
sides  of  the  TE  modules  increases.  The  same  remarks  relating  to 
minimum  and  higher  values  of  ATcoid-hot  may  be  repeated. 

Fig.  13  illustrates  the  power  trend  of  the  hydraulic  pump  vs 
speed  of  the  working  fluid  as  a  function  of  the  inner  diameter  as 
well  as  of  the  surface  area  of  the  solar  collector.  As  expected,  the 
power  of  the  hydraulic  pump  increases  with  increasing  of  speed 
of  the  working  fluid  (see  Eqs.  2,  3  and  11),  of  the  inner  diameter 
(see  Eqs.  1,  2  and  11)  as  well  as  of  the  surface  area  of  the  solar 
collector  (obviously  major  and  minor  losses  increase). 
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Fig.  5.  Net  electrical  energy  production  of  the  whole  system  in  the  case  of  collector  surface  equal  to  8  cm2  and  16  m2,  for  various  diameters  of  the  tube. 


Fig.  6.  Net  electrical  energy  production  of  the  whole  system  in  the  case  of 
D,  =  9.6  mm,  by  varying  the  collector  surface. 


Speed  of  working  fluid,  u  (m/s) 

Fig.  7.  Maximum  efficiency  of  the  whole  system  in  the  case  of  D,  =  9.6  mm,  by 
varying  the  collector  surface. 

In  order  to  assess  the  magnitude  of  the  economic  feasibility  of 
the  proposed  system  utilizing  the  Peltier  cells,  a  rough  estimate 
cost  per  unit  of  energy  produced  was  performed. 

The  cost  of  a  single  thermocouple  was  estimated  at  around 
0.22  €/thermocouple  considering  a  mean  current  commercial  cost 


Fig.  8.  System’s  maximum  electrical  power  in  the  case  of  D,  =  9.6  mm  varying  the 
collector  surface. 


Fig.  9.  Maximum  temperature  difference  of  working  fluid  in  the  case  of  D,  =  9.6 
mm,  by  varying  the  collector  surface. 


of  available  thermoelectric  modules.  The  cost  of  a  thermal  solar 
collector  was  estimated  at  around  800  €/m2  on  mean  market 
prices.  The  overall  cost  of  the  pipes  and  hydraulic  components 
was  estimated  at  around  500  €  on  mean  market.  Considering  the 
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Fig.  10.  The  trend  of  temperature  difference  between  two  sides  of  the  TE  modules, 
for  D,  =  9.6  mm  and  Ac  =  8  m2. 


Outside  air  temperature  (°C) 


0  10  20  30  40  50 

Outside  air  temperature  (°C) 


Outside  air  temperature  (°C) 

Fig.  11.  The  temperature  difference  between  two  sides  of  the  TE  modules  versus 
hourly  outside  air  temperature  (for  three  different  values  of  speed:  a.  u  =  2.3,  b. 
u  =  1.3,  c.  u  =  0.3),  for  D,  =  9.6  mm  and  Ac  =  8  m2. 

case  of  a  total  thermal  collector  surface  equal  to  16  m2,  a  pipe  with 
the  inner  diameter  of  14  mm,  therefore  7000  thermocouples  for 
10  m  of  length  of  thermoelectric  device,  a  working  fluid  speed  of 
0.8  m/s,  we  roughly  obtained  an  installation  cost  of  around 
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Hourly  horizontal  global  solar  radiation  (MJ/m2) 
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Fig.  12.  The  temperature  difference  between  two  sides  of  the  TE  modules  versus 
hourly  global  solar  radiation  (for  three  different  values  of  speed:  a.  u  =  2.3,  b.  u  =  1 .3, 
c.  u  =  0.3),  for  D,  =  9.6  mm  and  Ac  =  8  m2. 


Speed  of  working  fluid,  u  (m/s) 


Fig.  13.  Trend  of  hydraulic  pump  power. 

800  x  16  +  0.22  x  7  000  +  500  =  14  840.00  €  with  an  annual  net 
energy  production  of  197.27  kW  h/a.  Considering  a  system  lifetime 
of  25  years,  we  get  a  cost  of  3.01  €/l<W  h. 

This  cost  is  quite  expensive  compared  to  the  corresponding  cost 
of  the  photovoltaic  system.  Indeed,  considering  a  PV-system,  we 
have  a  mean  installation  cost  of  around  2000  €/l<W  p  and  the 
related  expected  energy  production  in  the  South  of  Italy  is  at 
around  1500  kWh/a.  Considering  a  system  lifetime  of  25  years, 
we  evaluated  a  cost  of  0.05€/kWh,  a  value  56  times  less  than 
the  proposed  system. 
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5.  Conclusions 

The  proposed  system  utilizes  the  thermal  difference  between  a 
lower  temperature  sink  and  a  higher  temperature  source.  The  sys¬ 
tem  is  focused  on  the  thermoelectric  generator,  basically  a  heat 
exchanger  composed  of  a  set  of  thermocouples.  The  main  parame¬ 
ters  of  the  system  are  the  thermal  power  and  the  temperature  of 
the  heat  source.  Considering  these  parameters  and  other  character¬ 
istic  variables,  simulating  software  was  developed  by  the  authors 
to  perform  several  tests. 

Results  show  that,  with  reference  to  the  annual  energy  produc¬ 
tion,  an  optimization  of  the  system  is  possible  choosing  a  suitable 
set  of  design  parameters,  basically:  an  inner  diameter,  a  solar 
collector  surface  and  a  working  fluid  speed.  As  expected,  the 
outcomes  also  point  out  that  the  size  of  the  inner  diameter  of  pro¬ 
posed  system  depends  on  amount  of  available  heat  flux  and  that 
the  efficiency  of  the  proposed  system  is  less  than  5%,  considering 
the  actual  material  thermocouples. 

The  complexity  and  the  poor  efficiency  of  the  proposed  system, 
as  well  as  the  high  cost  of  materials,  discourages  this  type  of  appli¬ 
cation  of  thermoelectric  generators. 

However,  the  results  obtained  suggest  further  investigation  on 
newer  thermoelement  materials  in  order  to  achieve  higher  effi¬ 
ciency  and  lower  costs. 
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